All-optical control of the spin state in the NV -center in diamond 
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We describe an all-optical scheme for spin manipulation in the ground-state triplet of the nega- 
tively charged nitrogen-vacancy (NV) center in diamond. Virtual optical excitation from the 3 A2 
ground state into the 3 E excited state allows for spin rotations by virtue of the spin-spin interaction 
in the two-fold orbitally degenerate excited state. We derive an effective Hamiltonian for optically 
induced spin-flip transitions within the ground state spin triplet due to off-resonant optical pumping. 
Furthermore, we investigate the spin qubit formed by the Zeeman sub-levels with spin projection 
ms = and ms = — 1 along the NV axis around the ground state level anticrossing with regard to 
full optical control of the electron spin. 



I. INTRODUCTION 

Nitrogen-vacancy (NV) centers in diamond have at- 
tracted much attention in research related to quantum 
computation^ due to their key advantages, such as high 
stability and long spin coherence timea^EI up to room 
temperature and beyonoP. The spin coherence time can 
be increased further by isotopic engineering^ since only 
the 13 C carbon atoms have non-zero nuclear spin, thus 
contributing to spin decoherence due to hyperfine cou- 
pling. Under resonant optical excitation the NV~ center 
exhibits a strong and highly stable zero phonon line at 
1.945 eV^ with an excited state lifetime of about 12 ns^. 
Electron spin resonance analysis of the center has shown 
that both ground state and excited state are spin triplets, 
which implies that there is an even number of active 
electrons involved. The ground state levels with spin 
projection TO5 = and mg = —1 along the NV axis 
become degenerate in a magnetic field of about 1025 G. 
Optical pumping causes a spin polarization of the ground 
state^"^ that can be attributed to a spin-orbit induced 
intersystem crossing with an intermediate sing let stattpl. 
When the zero field splitting is larger than the optical 
linewidth, repeated optical excitation leads to a spin se- 
lective steady state population in the lowest ms — level 
of the ground state, generating a non-Boltzmann steady 
state spin alignment and mixing of spin states^, so the 
spin of the ground state can be both initialized and read 
out optically^. 

The standard procedure for spin manipulation in 
the ground state triplet involves an oscillatory (radio- 
frequency) magnetic field that gives rise to electron spin 
resonance. In this paper, we describe an alternative 
method for full spin control without rf-fields, based en- 
tirely on optical transitions. All-optical spin manipula- 
tion of NV centers could allow for fast operations with 
high spatial resolution. In semiconductor quantum dots, 
picosecond optical control of single electron spins has 
been achieve d 13 ^ I Optically induced spin rotations in 
a single NV center in diamond have been demonstrated 
using off-resonant laser excitation^. This type of spin 
control relies on the optical Stark effect, i.e., the shift in 
energy levels induced by an applied optical field. Here, we 
describe an extension of this scheme which also allows for 




Figure 1. NV-center in diamond with equivalent sp 3 - 
hybridized dangling bonds (blue) Oi (i = 1,2,3) from the 
surrounding carbon atoms (black) and <tjv from the substitu- 
tional nitrogen atom (orange) overlapping within the vacancy 
(center). The NV-axis (long red arrow) defines the z-axis 
of the coordinate system and, due to the zero field splitting 
determines the quantization axis of the electron spin. The 
magnetic field B is applied along the NV-axis. 



transitions between the three ground state levels, similar 
to existing schemes for coherent population trapping^. 



II. MODEL 

To model optical spin rotations in an individual NV- 
center in diamond, we start from the commonly used 
description that fundamentally involves a total number 
of six e lectrons , but can be reduced to an effective two- 
electrori- 12 * 19 * 20 ! or, equivalently, two-hole model 21 . The 
four relevant single-electron orbitals ai, 02, e Xl and e y 
can be obtained by projecting the sp 3 -hybridized dan- 
gling bonds 01,2,3 of the carbon atoms and <jn of the 
nitrogen atom (see Fig. [T]) onto the irreducible represen- 
tations of the Cz v symmetry group of the NV center^. 
The electron configurations for the ground and first ex- 
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cited state are obtained as follows: In the ground state 
(e ) configuration, the lower-energy ai^ orbitals are com- 
pletely filled with two electrons each, while the e x and 
e y contain one electron each. The two-fold degenerate 
excited state configuration (ae) is obtained by promot- 
ing another electron from a\ to e x or e y . Due to the 
Coulomb interaction between the two electrons, the spin 
triplet lies lowest in energy and forms the ground state 
e 2 (T), transforming according to the representation A 2 
of Electric dipole transitions connect this triplet 

to the excited state spin triplet ea(T), transforming ac- 
cording to the E representation. The two-fold orbital 
and three-fold spin degeneracies give rise to a total of six 
states in ea(T), compared to three states in e 2 (T). The 
spin singlet states will not be of direct importance for 
our discussion, and are left out of our model. The entire 
state space for our model is thus nine-dimensional. 



where B denotes an external magnetic field aligned with 
the iVF-axis, g gs the Lande g-factor, and \x B the Bohr 
magneton. Around the ground-state level anticrossing 
(LAC), we can split the Zeeman energy into a term that 
compensates the ground state zero field splitting and an 
additional variation, g[i B B = D gs + g^i B 8B. The zero 
field splitting D gs = 2.88 GHz 22 is caused by the reduc- 
tion of the symmetry in spin space to C^ v due to the 
crystal field. The absence of orbital degeneracy in the 
ground state triplet implies that strain and spin-orbit in- 
teraction have very little effect on the ground state. 

Here, we have neglected the effect of the hyperfine cou- 
pling to the nuclear spins of the intrinsic nitrogen atom 
and surrounding 13 C atoms. If necessary, the nuclear 
spin state could be prepared optically. 



A. Ground state 

The Hamiltonian of the ground state spin triplet in the 
basis { 3 A 2 _, 3 A2o, 3 A 2+ } = {| - 1), |0),| + 1)} is 

f D gs - g gs ^ B B 

//... = I I , (1) 

D gs + g gs /j, B B, 



B. Excited state 

At low temperatures, the excited state fine structure 
can be understood to a large extent from strain and 
spin-spin interactions. In the basis of spin-orbit states 
with full symmetry, described by the Cz v double group 
including spirit, {Ai,A 2 ,E x ,E Y ,E 1 ,E 2 } the excited- 
state Hamiltonian matrix is 
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where l z is the axial spin-orbit splitting, and D cs = 
— jD zz and A = ^D x 2_ y 2 are the well-known spin-spin 
interactions^. Experimentally, it was found that l z = 5.3 
GHz, D CB = 1.42 GHz and A = 1.55 GHsP. The Lande 
factors of ground and excited state were found to be 
equapgt!!] ; ^ gs ~ 5os ~ 2.01 = g. The energy gap E g is 
defined as the difference between the E\^ 2 excited states 
and the rag = ±1 ground states at B = 0. The de- 
pendence of the ground- and excited state levels on B is 
shown in Fig. [2] 

Since electric dipole transitions are spin-conserving, 
our all-optical spin control scheme requires a spin non- 
conserving mechanism in the excited state. The longitu- 
dinal spin-orbit interaction term l z only leads to an addi- 
tional energy splitting between states with different spin 
projections and cannot flip the spin. It was speculated 
that the transv ersal p art of the spin-orbit interaction can 
lead to spin flip d 23 * 24 ! , but it has recently turned out that 
it can only connect orbital states belonging to different 
irreducible representations^. However, the transversal 
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component A" ~ 0.2 GHz of the spin-spin interaction al- 
lows for the non-spin-conserving transitions between the 
E and E' states, explaining the experimentally observed 
transitions^. 

The components S x and S y of the non-axial strain 
can be written in polar coordinates, 5\ = 5 X + S y and 
2/3 = arctan^j,/^). Here, j3 was defined such that it 
corresponds to the angle between the symmetry axis of 
strain eigenstates e' x (S Xl S y ) and the symmetry axis of 
the unperturbed e^-orbital. 



C. Electric dipole transitions 



We assume the system to be optically driven with a ra- 
diation field at fixed frequency u> near E g , therefore it is 
convenient to describe the excited states in a corotating 
frame, while keeping the ground states fixed. We then 
work in the rotating wave approximation where counter- 
rotating terms with frequency E g +uj are neglected. This 
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and (i = x, y) 

e, = MJMM^ (7 ) 

with the reduced matrix element of the position operator 
defined as 

(e| \f E \ \a) ={e x \x \a) = (e y \ y \a) . (8) 

We are interested in linearly polarized optical fields, 
where e is real. The magnitude of the dipole matrix 
elements can be estimated from the observed Rabi os- 
cillation periocP3 and the linear Stark shiftP^, typically 
e re 1 GHz. 



III. EFFECTIVE SPIN HAMILTONIAN 
A. Schrieffer- Wolff transformation 



Figure 2. Illustration of the optical pumping process with 
frequency to from m s = ground state spin level below the 
excited state. The degenerate E x , v levels are coupled to the 
-Ei, 2 levels by the transversal spin-spin-interaction A". 

is justified as long as 5u) = E g — uj <C E g . Optical 
transitions between the ground and excited state are de- 
scribed with the electric dipole operator for two electrons, 

tfS=2T$®l + l®H$> where 

H^l = eE • f = e|E| (xcosa + ysina) (3) 

is the single-particle electric dipole operator, where f = 
(x, y, z) denotes the electron position operator. Here, 
we assumed the incident light to be linearly polarized 
perpendicular to the NV-axis (which defines the z-axis 
of the coordinate system) with polarization angle a, with 
a = for polarization parallel to the symmetry axis of 
the e x -orbital. Taking matrix elements with the ground 
and excited state basis states yields the Hamiltonian 

where the detuning Suj is defined with respect to the 
lowest-lying excited state energy levels Ei t 2 (neglecting 
strain and spin mixing A"), and Auj = 8lu + l z — 
The transition matrix is given as 

(ie + — it- — 

-1e y 2e x (5) 
— ie_ — ie_ ie + J 

where, for linear polarization of the excitation field, 

e± =e x ±ie y = ee ±la , (6) 



Since the energy levels of the ground state H gs and 
the excited state H cs are widely seperated by E g = 1.945 
eV sa 470 THz and coupled by small perturbations 
e -c \6(jj\ = \E g — ui\, we can use a Schrieffer- Wolff 
transformation's of the Hamiltonian (block-matrix) as 
a valid approach to determine the effective dynamics of 
the driven system up to second order in the perturbation. 
The Schrieffer-Wolff transformation is defined as follows, 

H = e s He- s = H+[S,H} + ^ [S, [5, H]]]+0 (S 3 ) , (9) 
with the anti-hermitian transformation matrix 

»-(_ .,;). <m 

The aim of the transformation is to remove the coupling 
V in first order, which can be achieved if [S, H ] = — V. 
In terms of the submatrices for the two seperated systems 
H gs and H es and their coupling V this condition reduces 
to 

sH cs - H gs s = -v, (11) 

which also implies that in a perturbation series in v, the 
matrix s will be first order, s ~ 0(v). This particular 
choice of transformation secures that the first order terms 
in Eq. ^ cancel and we are left with an effective ground- 
state Hamiltonian 

H = H gs + i (sWAs) +0(v 3 ) . (12) 

Note that in the absence of strain l z — D cs — D gs s» 1 GHz 
is the seperation between the two closest-lying energy 
levels E XtV and E + = Ei + E%i and thus for resonant 
excitation between those two levels, the optical driving 
field strength has to be much smaller than 500 MHz. 
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Figure 3. The state |s) of the ground state spin qubit in the 
ms = 0,-1 subspace can be represented as a vector on the 
Bloch sphere where the two poles represent the two Zeeman- 
split eigenstates \m s = 0) and \m s — —1). The vector b de- 
notes the effective magnetic field acting on this pseudo-spin 
1/2, with spherical angles 9 and </>. 



B. Rotation axis 

Wc focus on the transition between the ms — — 1 and 
ms — ground state spin levels. This two-level system 
can be split off from the ms = +1 state^, in the regime 
g^B\SB\ < 2Z? gs . In this case the dynamics is described 
by the 2x2 Hamiltonian 



H = 



Hll #12 
if 21 H22 



= 6 H + b • er 



(13) 



where the effective (pseudo-) magnetic field has compo- 
nents, 



b x = i (J/12 U 2 i ) = b_ nisi.. 



(14) 



b y = -(H 12 -H 21 ) =b ± smq>, (15) 



K = 



(16) 



and bo — -Hn + Hti- The dynamics within this two- 
dimensional subspace can be visualized using a Bloch 
sphere picture, as shown in Fig. [3| Optical driving re- 
sults in a rotation of the state vector about an axis 
b with polar angle (axial orientation) 9 and azimuthal 
angle (non-axial orientation) <f>. The transversal part 



bl + bl 



b sin 9 of the effective field pro- 
vides for effective spin-hip transitions while b z accounts 
for the effective (AC Stark) splitting between the two 
levels (around the LAC). The term cx 1 in Eq. (13) can 
be omitted since it merely leads to a global phase. The 
precession frequency is given by 



b = Jb\ + b\ = y/l + cot 2 9\b ± \. 



(17) 



Spin flips can be implemented as rotation about an axis 
within the equatorial plane of the Bloch sphere (Fig. [3]) 
which corresponds to 9 = f and thus b = bx ■ 



IV. RESULTS 



Unstrained NV center 



In the case of vanishing strain (5± — 0) we obtain a 
simple analytical result for the transversal component of 
the qubit rotation axis, with magnitude, 



b x = A"e 2 



1 



Sw (Suj + D cs + l z + -Dgs + g^B&B) 



(18) 



(5u - gfJ-B^B) (5uj - D cs + l z + D gs ) 



0(A 



l/2\ 



which is proportional to the intensity e 2 of the optical 
driving field and the transversal spin-spin coupling A" 
in the excited state. The azimuthal angle of the rotation 
axis is determined by the optical polarization angle a, 



-2a 



(19) 



where the factor of 2 reflects the double group character 
of spin representation. The polar angle 9 of the rota- 
tion axis is independent of a and for small g^isSB even 
independent of the driving field strength e. The resid- 
ual Zeeman splitting is limited by the hyperfine LAC of 
about 2 MHzS3, and therefore for an optical coupling 
e > 100 MHz we can always find pairs of parameters 
(SB, Sw) that fullfill the condition 9 = tt/2. 

In the limit of large detuning, i.e., when Suj dominates 
all other energies in the denominators of Eq. ( |18[ ) , we can 
approximate the transverse component of the effective 
field as 



fe ± ~ 2A" 



(20) 



B. Effect of strain 



We now include the effect of strain in the diamond 
crystal into our discussion. For moderate strain S XtV <C 
2l z m 10 GHz, the levels of the excited state are 
largely seperated from the E\ t i levels, and thus the 
strain-induced mixing of states and states can 
be neglected in lowest order. The main effect of moder- 



ate strain is thus a shift of the resonances in Eq. ( 18 ) by 

±Sx = ± ^Js 2 + <5 2 , lifting the degeneracy of E x and E y 

levels. Though strain does not directly mix states with 
different spin projections, this shift reduces the energetic 
seperation between coupled E and E' levels and therefore 
strongly enhances the efficiency of spin-flip-transitions. 
In Fig. |4j we plot suitable pairs of parameter values for 
the detuning Suj (near resonant driving) and Zeeman 
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Figure 4. Pairs of required magnetic field variation SB and 
detuning 5uu at fixed dipole coupling strength e = 500 MHz 
for different values of transversal strain S± — 1,2,3,4 and 5 
GHz (from top to bottom) to fullfill the condition 6 = 5 for 
precession around a rotation axis within the equatorial plane 
of the Bloch sphere for a = ft = 0. 
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Figure 5. Precession frequency b± for different values of 
transversal strain 5± = 1,2,3,4 and 5 GHz (from left to 
right) for polarization parallel (a = ft) and perpendicular 
(q — P = 5) to strain with dipole coupling e = 500 MHz. 



Figure 6. Pairs of detuning &ui and transversal strain S± that 
match the condition 9 = \ for precession around a rotation 
axis within the equatorial plane of the Bloch sphere for differ- 
ent dipole coupling e = 100, 200, 400, 600, 800 and 1000 MHz 
(from bottom to top) for a = ft — 0. 
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Figure 7. Precession frequency fox around in-plane rotation 
axis at 10 GHz detuning for optical coupling e = 500 MHz 
and polarization in e x (a = 0) and e y (a — ?) direction 
for different values of transversal strain S± = 1, 2, 3, 4, 5 GHz 
(increasing amplitude) over strain direction angle ft. 
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Figure 8. Non axial orientation angle </> of precession axis for 
polarization parallel (a = 0) and perpendicular (a = f ) to 
e^-orbitals symmetry axis in respect to strain orientation for 
different strengths of transversal strain. 



splitting gfiB^B with varying strain 5j_ that fulfill the 
condition 9 = \ for an in-plane rotation axis. This de- 
fines an implicit function SB (8w, 8±) given a and /3, e.g. 
in the case of 9 = tv/2 for a = j3 = 0. To find the strength 
of spin-flip transitions, we substitute SB into the preces- 
sion frequency for an in-plane rotation axis, b± and plot 
it in Fig. [5] as a function of the optical frequency detun- 
ing Slu for different values of transversal strain. We find 
that the precession frequency b± indeed increases with 
strain. Varying e numerically shows that the precession 
frequency is still proportional to the intensity of the op- 
tical driving field 6j_|e=f oc e 2 . Note that the perturba- 
tive approach breaks down as detuning approaches strain 
(divergence of bj_ in Fig. [5]), restricted by the validity 
condition Slu <C <5j_ for the Schrieffer- Wolff transforma- 
tion in Eq. Q. We also investigate the dependence on 
the direction of strain and polarization (see Fig. [7f. Ex- 
pectedly we get the highest efficiency for collinear strain 
and polarization and a minimal efficiency for perpendic- 
ular relative orientation with an overall sinusoidal form 
of twofold symmetry. Changing the optical polarization 
angle a only leads to a uniform and continuous shift of 
this function (this has been checked for a variety of dif- 
ferent values, but for simplicity we only show it for a = 
and a = Thus, for weak strain, the resulting effective 
field only depends on the relative angle a — j3 between the 
strain and polarization angles. However, as the transver- 
sal strain S± increases beyond about 10 and 20 GHz, we 
start observing a modulation of the field with higher har- 
monics of a — j3. Numerical evaluation also reveals that 
the azimuthal angle (for in-plane orientation of the pre- 
cession axis) is generally independent of optical coupling 
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Figure 9. Non axial orientation angle (j> of precession axis for 
different polarization angles a in respect to strain orientation 
at transversal strain 5± = 5 GHz. 

strength and magnetic field, 

4> = 4>(a,P,6j_,5u). (21) 

From Fig. [8] we see that for e x /j,-polarized light (i.e. 
a = 0,5), the angle <fr is well approximated by <p = 

- s r4' 5uj) (5±, Slu) sin 2/3 (at low strain), with an am- 
plitude proportional to the intensity of the strain, 

= 6 ± f (Su) (22) 

In Fig. [9j we show that the sinusoidal shape of as a 
function of a is slightly distorted for polarization angles 
a =/= ri^, n G N; we also find that this distortion grows 
with increasing strain. 

C. High strain limit 

In the high strain regime, the electronic states of 
the NV-center are energetically split into two orbital 
branches with largely seperated energies E x and E y cor- 
responding to a specific choice of coordinate axes, that 
fixes the orientation angle (3. In this limit b± and b z (and 
thus the polar angle 0) become independent of the orien- 
tation angles of both strain and polarization. Expanding 
in 5 J 1 yields 

b ± ~ 2A"|-, (23) 

and 

4 = 2(a-2/3)-^ + o(±-Y (24) 

for the transversal part of the pseudo-field, and 

^-^ + 2(^-^1 + 0^), (25) 

for the longitudinal part. The higher orders contain 
higher harmonics, such as terms cx A cos 2 (a — 3/3) for 

b x and oc A sin 2 (a — 3/3) for b y in O^jt^j, indicating 

strain induced third order transitions mediated by the 
Ai.2 levels. 
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CONCLUSION 



We have shown that the effective precession axis and 
frequency of the ground state spin of the iW~-center 
can be fully controlled by off-resonant optical excitation, 
by adjusting the frequency detuning Suj and linear po- 
larization angle a of the optical driving field for a given 
intensity e (optical dipole coupling) and magnetic field 
B. The orientation of the precession axis is determined 
by two angles 9 and (f>, where the first is depending on all 
parameters (including strain and polarization) and the 
latter is independent of magnetic field and optical cou- 
pling strength and basically controlled by polarization 
and strain. The strain effects can be compensated by ex- 
ternal bias voltage^. Since any unitary qubit operation 
(rotation around axis n by angle 7) can be composed by 
successive rotations around two orthogonal axes on the 
Bloch sphere, a complete set of single-qubit operations 
can be generated optically in this way. From a purely ge- 
ometric point of view, spin rotation about an axis within 
the equatorial plane of the Bloch sphere (where b z — 0) 
is most effective for flipping the spin, although any axis 
other than the z-axis would do (the smaller the polar an- 



gle 9, the more pulses are be required). A full spin- flip 
is obtained by a b±r — 7r-rotation around an axis within 
the equatorial (x, y)-plane of the Bloch sphere, provid- 
ing an estimate for the gate switching (optical pumping) 
time r in the limit of large detuning and weak strain, 



7r Suj 2 



2 A"e 2 ' 



(26) 



The switching time for spin-flip transitions is limited by 
the spin mixing term A", since the (above) condition im- 
plies via the off-resonant condition e <C Sui the following 
lower limit for the spin- flip, 



r > — w 10 ns 
~ A" 



(27) 



where for that latter estimate we assumed A" = 0.2 GHz. 
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